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ABSTRACT 


Concurrent  water  samples  and  optical  groflles  gere  taken  at  twelve  stations  ln‘the 
Western  Sargasso  Sea,, at  approximately  34°N  and  70  W  (Blowatt  Mooring)  during  August  1987. 
v__»>0ptTcaI "measuremSnCs  included  up-  and  downwelling  scalar  and  vector  Irradiance.  In  situ 
relative  fluorescence  and  beam  transmission  at  660  nm  revealed  a  particle  layer  the  center 
of  which  varied  between  70  emd  110  m  and  had  a  breadth  of  20  to  40  m.  Water  samples  were 
collected  within  the  particle  layer  and  near  surface  and  analyzed  for  particle 
characterization  and  volume  absorption. 


Absorption  coefficients  for  the  near  surface  and  particle  layer  were  calculated  using 
-threa_ techniques.)  Two  techniques  employed  measurements  of  apparent  optical  properties. 

>0ne  method  utilized  scalar  and  vector  irradlances  and  calculated  absorption  as  the  product 
of  the  average  cosine  and  the  diffuse  attenuation  coefficient  for  net  irradiance  at 
wavelengths  of  441  and  520  nm.  The  second  method  used  inversion  of  the  irradiance  quartet 
(up-  and  downwelling  vector  and  scalar  irradlances)  and  also  calculation  of  absorption  from 
only  up-and  downwelling  vector  irradlances.  In  the  third  method,  absorption  coefficients 
were  calculated  as  the  sum  of  particulate,  dissolved  substances,  and  molecular  water 
absorption.  Some  large  discrepancies  were  found  between  absorption  from  apparent 
properties  and  those  calculated  from  the  component  sum.  Possible  sources  of  error  include 
overestimation  of  dissolved  substance  or  particulate  absorption,  inaccurate  ''^re  water 
absorption,  and . inf luence  of  internal  radiance  sources,  such  as  Raman  scattering,  on 
apparent  optical  properties.  /  '5  /’//  _ 

Particles  were  characterized  by  three  (different  techniques:  relative  particle 
concentration  by  in  situ  transmissometry  at  660  nm;  particle  frequency  and  volume 
distributions  by  resistive  pulse  technique  (’Coulter  counts”);  and  particle  frequency  and 
size/shape  distribution  vs.  elemental  composition  by  SAX  (Scanning  electron  microscopy  with 
Automated  image  analysis  and  X'l'sy  spectroscopy).  Total  particle  number  decreased  with 
depth  and  was  largely  due  to  decreasing  concentration  of  organic  particles.  There  was  a 
small  size  increase  in  siliceous  particles  near  depths  of  the  fluorescence  amxima.  In 
contrast,  bacicseatter  values  derived  from  apparent  optical  properties  increased  by  nearly 
two-fold  from  the  surface  to  depth  of  maximum  fluorescence. 


1.  INTRODUCTION 


Absorption  and  scattering  are  fundamental  processes  governing  light  attenuation  in 
natural  waters.  Preisendorfer'*  has  distinguished  between  inherent  optical  properties,  such 
as  the  volume  absorption  and  scattering  coefficients  which  are  functions  of  the  medium  and 
are  Independent  of  the  directional  nature  of  the  light  field,  from  apparent  optical 
properties,  such  as  the  diffuse  attenuation  coefficient  (K)  and  irradlances  which  are 
dependent  on  the  light  field.  There  are  no  accepted  awthods  to  independently  measure 
absorption  and  scattering  coefficients  in  situ;  however,  absorption  and  scattering  can  be 
estimated  frewn  laboratory  measurements  of  discrete  water  samples  or  they  can  be  derived 
through  inversion  of  radiative  transfer  models  which  use  apparent  optical  properties  such 
as  measured  up-  and  downwelling  irradlances.  A  goal  of  the  Naval  Ocean  Research  and 
Development  Activity  fNORDA)  optical  oceanography  program  is  deriving  inherent  optical 
properties  for  various  ocean  environments,  then  checicing  the  volume  or  "bullc”  optical 
properties  for  internal  consistency  using  Monte-Carlo  simulation. 

Even  though  various  techniques  to  estimate  the  inherent  optical  properties  have  existed 
for  some  time,  few  investigators  have  attempted  to  compare  laboratory  measurements  of  bulk 
properties  with  inherent  properties  derived  from  theoretical  zK>dels  that  use  apparerit 
properties  as  inputs.  Ne  present  optical  and  particle  data  from  which  the  absorption 
coefficient  is  estimated  for  blue  and  green  light  from  three  methods:  laboratory 
Masureamnt  of  discrete  water  samples,  inversion  of  in  situ  irradiance  measurements,  and 
solution  of  the  radiative  transfer  equation  for  net  irradiance. 
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2.  THEORY 


2.1»  Absorption  via  Stavn 
2 

Stavn  expanded  on  the  classic  work  of  Gershun  and  Schellenbecger  to  derive  what  he 
termed  the  three-parameter  model  for  a  light  field.  In  his  work^  the  absorption 
coefficient  for  a  depth  Interval  is  calculable  from  measurements  of  up-  (E  )  and 
downwelllng  (E.)  vector  and  scalar  (E  )  irradiances.  Whereas  other  treatments  have  largely 
ignored  changes  in  the  average  cosine°(mean  light  path)  that  may  occur  over  a  depth 
interval,  Stavn  showed  that  this  may  lead  to  an  overestimate  in  absorption.  The  absorption 
coefficient  (a)  for  a  depth  interval  to  z^  is  given  for  the  interval  by 

a  ■  (In  E^(Z2)  -  In  E^Cz^))/  A2  (1) 

where  E.  *  E .  -  E  »  net  irradiance 
Z  Cl  u 

AZ  -  dz/u  -  jQ^l  dz/u 

u  *  average  cosine  given  by  (E^/E^) 

The  mean  light  path  (Z  above)  retains  the  information  regarding  variation  in  the  average 
cosine  or  changes  in  light  path  occurring  within  the  depth  interval  over  which  absorption 
is  being  calculated.  Equation  (1)  can  be  written  as 

a  -  u  Kg  (2) 

from  the  definition  of  the  diffuse  attenuation  coefficient  for  net  irradiance  (K  )  and  with 
u  averaged  over  the  s^m|  space  and  time  interval  from  which  K  is  calculated.  T^e  Stavn 
three-parameter  model  '  is  limited  in  that,  while  the  absorption  coefficient  can  be 
extracted,  the  scattering  coefficient  cannot. 

2.2.  Absorption  via  Preisendorfer/Mobley 

By  contrast,  both  absorption  and  backscattering  coefficients  can  be  obtained  from 
inversion  of  the  irradiance  quartet  (up-  and  downwelllng  scalar  and  vector  irradiances). 
Preisendorfer  and  Mobley  have  suggested  that,  for  most  natural  waters  the  upwelling 
(E  /E  )  and  downwelllng  (E./E  .)  average  cosines  are  approximately  constant  with  values  of 
l/1.5x  and  1/1.33,  so  that  the  volume  absorption  a  and  volume  backscatter  b.  can  be 
calculated  simply  from  the  up-  and  downwelllng  vector  irradiances.  “ 

2.3.  Absorption  by  components 

Because  absorption  is  an  inherent  property,  it  can  be  partitioned  into  its  components  so 
that  the  contribution  from  water  a  ,  particulates  a  ,  and  dissolved  substances  a  can  be 
evaluated.”''.  Thus,  P  y 


Inasmuch  as  K,  the  diffuse  attenuation  coefficient  iSgan  apparent  optical  property,  the 
same  type  of  breakdown  into  components  does  not  apply  . 

The  presence  of  internal  radiance  sources  such  as,  particulate  and  dissolved  material 
fluorescence,  Raman  scattering,  or  bioluminescence  can  significantly  influence  the  values 
of  absorption  derived  from  the  apparent  optical  properties,  but  in  general  they  will  have  a 
small  impact  on  the  absorption  calculated  independently  from  the  components.  Therefore,  a 
comparison  of  absorption  values  from  components  with  those  derived  from  apparent  optical 
properties  has  potential  for  quantifying  the  combined  effects  of  bioluminescence, 
fluorescence,  and  Raman  scattering. 


3.  METHODS 

3.1.  Sampling  area  and  mode 

Optical  and  particle  data  were  collected  at  one  station  along  the  north  wall  of  the  Gulf 
Stream  and  at  twelve  stations  in  the  general  vicinity  of  the  Biowatt  mooring  at  "34°  N  and 
72°  W.  We  report  here  preliminary  results  for  two  stations  near  the  Biowatt  mooring. 
Stations  8  and  11.  Deployment  at  Station  8  (21  August)  was  a  typical  vertical  cast  whereas 
Station  11  (22  August)  consisted  of  3-5  minute  time  series  at  depths  of  10,  30,  65,  75,  and 
100  m. 
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tlcal  and  hydrographic  mcasureaents 


Optical  data  and  water  sanples  were  collected  with  a  Particle/Optical  Sampling  System 
(POSSY),  consisting  of  in  situ  instrumentation  attached  to  a  General  Oceanics  (Miami,  PI.) 
rosette  with  eight  liter  Niskin  water  samplers.  The  optical  suite  consisted  of  a 
Biospherical  Instruments,  Inc.  (San  Diego,  Ca.)  MER-1048  spectroradiometer  which  measured 
down-  (13  wavelengths)  and  upwelling  (8  wavelengths)  vector  irradiance,  upwelling  radiance 
(8  wavelengths),  and  up-  and  downwelling  scalar  irradiance  (4  wavelengths  each).  Due  to 
the  loss  of  a  downwelling  scalar  sensor,  the  complete  optical  quartet  was  only  available 
for  441,  488,  and  520  nm. 

Beam  transmission  was  measured  with  a  25  cm  path  length  Sea  Tech,  Inc.  (Corvallis,  Or.) 
transmlssometer,  and  relative  fluorescence  was  measured  with  a  Seamartec,  Inc.  (Seattle, 
Wa.)  in  situ  fluorosMter.  Outputs  from  both  instruments  were  multiplexed  into  the  PSK 
signal  sent  to  the  surface  by  the  HER.  Temperature  and  conductivity  were  measured  using 
Sea-Bird  Electronics  Inc.  (Bellevlew,  Wa.)  modular  sensors  whose  outputs  were  also 
multiplexed  into  the  PSK  output. 

Water  samples  were  collected  on  upcasts  or  during  time  series  for  particle  analysis 
(elemental  composition,  sise,  concentration,  volume),  and  dissolved  material  and 
particulate  absorption  spectra.  Dissolved  absorption  was  measured  from  320  to  760  nm  on 
the  filtrate  passing  through  a  0.2um  nucleopore  filter  using  a  Kontron  Uvikon  Model  860 
(Everett,  Ma.)  spectrophotosMter.  Because  dissolved  substance  absorption  was  so  small,  the 
natural  log  of  absorbances  were  plotted  against  wavelength  from  320  to  420  nm,  and  the 
exponential  rate  of  decrease  was  calculated.  This  was  then  used  to  extrapolate  dissolved 
absorption  at  441  and  520  nm.  Particulate  absorption .was  determined  using  GP/C  filters 
following  procedures  put  forth  by  Mitchell  and  Kiefer  .  Particle  concentrations  and  size 
distributions  were  measured  by  two  techniques.  The  first  was  the  standard  "Coulter"  or 
resistive  pulse  technique  using  a  Particle  Data  Inc.  (Elmhurst,  II.)  Model  BO-XY  counter. 
The  second  was  an  automated  microscopy  technique  that  uses  a  scanning  electron  microscope, 
image  analyzer,  and  X-ray  spectrometer  combined  with  heavy  metal  staining  to  yield 
additional  information  on  particle  composition. 

3.3.  Irradiance  data  processing 

The  high  temporal  resolution  (5  Hz)  of  the  optical  data  allowed  for  the  resolving  of 
sources  of  contamination.  A  notch  filter  was  applied  to  the  raw  data  to  remove  irradiance 
fluctuations  that  were  linked  to  ship  motion  and  tilt-roll  of  the  instrument  package.  A 
technique  similar  to  that  used  by  Smith  and  Baker^  was  employed  to  smooth  out  effects  of 
wavefocusing  on  calculated  values  of  K  .  Because  its  variance  was  much  smaller  than  that 
of  K  ,  the  average  cosine  was  calculated  from  the  arthiieetic  mean  over  the  same  interval 
as  for  K.  To  account  for  temporal  changes  in  incident  light,  irradiances  were  normalized 
to  deck  cell  measurezMnts. 

For  the  Preisendorfer-Mobley  inverse  model,  irradiance  data  was  binned  within  +0.1  m  of 
each  integer  depth  for  profile  Station  8.  For  Station  11,  irradiances  at  each  integer 
depth  were  calculated  from  the  time  averages  at  each  of  the  five  depths,  normalized  to  the 
appropriate  deck  cell  (to  remove  atmospheric  variability).  These  values  were  then  used  to 
calculate  the  average  irradiances  at  intervening  depths  that  were  used  as  input  into  the 
Preisendorfer-Mobley  model.  The  application  of  splining  techniques  were  unsatisfactory  due 
to  the  limited  number  of  data  points. 

4.  RESULTS  AMD  DISCUSSION 


4.1.  Temperature,  fluorometrv,  transmissostetry 

Hydrographically,  the  mixed  layer  was  isothermic  at  "23.2  C  frosi  the  surface_to  35  m, 
from  which  depth  the  temperature  declined  at  a  steady  and  slow  rate  of  0.06  C  m~  to  at 
least  200  m.  The  beam  transmissosMtry  profile  for  Station  11  varied  +0.2%  about  88.3 
between  the  surface  and  200  m.  Relative  fluorescence  was  very  low,  vTrtually  at  zero  for 
most  of  the  water  column  except  between  150  to  165  m  and  between  65  and  110  m  (SEN 
examination  revealed  the  presence  of  chain-forming  diatoms  only  at  the  75  m  sample  depth). 
In  each  of  these  two  layers,  relative  fluorescence  rose  only  to  about  5%  on  the  most 
sensitive  scale.  The  beam  transmission  did  not  show  any  structure  at  the  more  shallow 
fluorescence  peak,  but  it  did  exhibit  a  small  decrease  at  the  deeper  maximum.  At  Station 
8,  the  in  situ  fluorescence  layer  was  between  80  and  120  m  peaking  at  "lOO  m._|At  either 
station,  extracted  chlorophyll  were  low  at  all  depths  rising  to  only  0.1  ug  1  in  the 
maximum. 
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4.2.  Absorption  coefficients 


Absorption  Measurements  are  calculated  for  two  wavelengths,  441  and  520  nm,  to  yield 
inforaation  on  the  spectral  dependency  of  bulk  absorption.  The  absorption  by  molecular 
water  a  for  441  nm  has  been  reported  to  be  0.014  ffl~  versus  0.048  m~  for  520  nm^  . 


Table  1. 

Bulk  absorpti 

on  of  water. 

dissolved  substances,  and  particulates  (m~^  ) 

Station 

8  441 

nm 

Absorption 

(m“^) 

depth 

(m) 

«w 

“p 

a 

a’ 

a" 

10 

0.014 

0.011 

0.013 

0.038 

0.012 

0.014 

0.014 

34 

0.014 

0.025 

0.013 

0,052 

0.017 

0.018 

0.018 

100 

0.014 

0.029 

0.039 

0.082 

0.053 

0.055 

0.051 

Station 

8  520 

nm 

depth 

(m) 

10 

0.048 

0.003 

0.004 

0.055 

0.042 

0.047 

0.040 

34 

0.048 

0.006 

0.003 

0.057 

0.039 

0.041 

0.035 

100 

0.048 

0.007 

0.008 

0.063 

0.041 

0.043 

0.052 

Station 

11  441 

nm 

depth 

(m) 

11 

0.014 

0.016 

0.012 

0.042 

0.014 

0.016 

31 

0.014 

0.009 

0.014 

0.037 

0.018 

0.016 

62 

0.014 

0.011 

0.020 

0.045 

0.029 

0.029 

80 

0.014 

0.015 

0.036 

0.065 

0.039 

0.044 

107 

0.014 

0.009 

0.015 

0.038 

0.038 

0.047 

Station 

11  520 

nm 

depth 

(m) 

11 

0.048 

0.004 

0.005 

0.057 

0.041 

0.043 

31 

0.048 

0.002 

0.004 

0.054 

0.040 

0.040 

62 

0.048 

0.002 

0.005 

0.055 

0.042 

0.042 

80 

0.048 

0.003 

0.011 

0.062 

0.044 

0.042 

107 

0.048 

0.002 

0.004 

0.054 

0.044 

0.043 

a'  »  aKsorp)flon  8erived  from  Preisendorfer/Mobley  using  calculated  distribution 
factors. 

a*  «  absorption  derived  from  Preisendorfer/Mobley  assuming  distribution  factors  of 
2.67  and  1.33  for  up-  and  downwelling  radiance  respectively. 
uK  »  absorption  calculated  using  Stavn's  three-parameter  method. 


4.3.  Absorption  by  dissolved  substances 

Stations  8  and  11  differ  markedly  in  the  absorption  attributed  to  dissolved  substances 
(Table  1);  however,  the  shape  of  the  absorption  spectra  were  quite  similar.  The  slope  of 
the  log-linear  relationship  between  In  of  absorption  and  wavelength  over  the  Interval 
320-420  nm  was  0.015-0.020  with  a  mean  of  0.017  at, both  stations.  This  range  for  the  rate 
of  extinction  is  the  same  as  that  found  by  others^'^'^  .  Using  absorption  at  380  nm  as  a 
reference  (a  (380)),  station  8  had  absorption  values  that  ranged  from  0.03  to  0.08  m  and 
station  11  had  values  less  than  0.05  m  .  The  a  (380)  at  other  stations  ranged  from  near 
detection  limits  to  "0.05  m~  .  .These  values  are^generally  in  the  range, found  by  other 
studies. for  oligotrophic  oceans^  (they  range  from  near  zero  to  0.05  m”  with  a  mean  of 
0.03  m"^). 

The  absorption  ascribed  to  dissolved  substances  showed  a  strong  relationship  to  particle 
absorption  and  in  situ  fluorescence.  Presumably  this  represents  a  non-conservative 
fraction  that  is  a  decay  product  of  phytoplankton'^  .  At  both  stations,  a  increased  by 
nearly  a  factor  of  2  near  depths  of  higher  relative  fluorescence.  ^ 

4.4.  Particulate  absorption 

The  magnitude  of  particulate  absorption  was  the  same  at  Stations  8  and  11.  Absorption 
ascribed  to  particulates  was  approximately  0.013  m~  near  the  surface  and  increased  nearly 
three-fold  in  samples  collected  at  or  near  the  fluorescence  maximum. 
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4.5«  Bulk  absorption  coeCflclents 
Station  8 

Station  8  was  a  vertical  profile,  and  exhibited  high  a  .  Irradiance  data  were  in  one  m 
increaents.  There  was  excellent  agreeaent  between  absorption  coefficients  derived  froa  the 
three-paraaeter  aodel  of  Stavn  and  the  Preisendorfer/Nobley  inverse  aodel,  which  assumed 
distribution  factors  (the  reciprocal  of  average  cosines)  of  1.33  and  2.67  for  down-  and 
upwelling  light  respectively  (Figure  1).  At  441  na  the  absorption  coefficients  derived 
froa  the_aodels  using  apparent  optical  properties  rise  froa  near  that  of  *pure  water* 

(0.014  a  '^)  at  the  surface  to  0.047-0.050  b~  at  the  fluorescence  maxima.  In  clear  waters, 
where  multiple  scattering  is  negligible,  it  was  anticipated  that  the  Stavn  and 
Preisendorfer/Nobley  models  would  yield  similar  results,  because  both  are  highly  dependent 
on  the  derivative  of  downwelling  vector  irradiance.  Including  the  scalar  Irradiances  and 
calculating  the  distribution  factors  (or  average  cosines)  rather  than  using  the  constants 
of  2.67  and  1.33  did  little  to  change  the  absorption  estimate  for  441  na  (see  Figure  2)  in 
the  Preisendorfer/Nobley  aodel. 

The  results  at  520  na  were  similar  in  that  the  Stavn  and  Preisendorfer/Nobley  models 
yielded  similar  absorption  estiatatesi  however,  this  was  only  true  if  the  distribution 
factors  were  calculated  from  the  irradiance  guartet  instead  of  using  1.33  and  2.67.  The 
Preisendorfer/Nobley  estimates  of  absorption  (using  the  calculated  distribution  factors) 
rose  froa  *‘0.040  above  50  a  to  >0.045  below  50  m;  whereas  the  absorption  coefficient  from 
the  three-parameter  aodel  was  *’0.046-0.050  a~  froa  10  to  70  a  then  declined  to  0.040  by 
100  m.  if  1.33  and_2.67  %rare  used  for  the  distribution  factors,  absorption  values  ranged 
froa  0.03  to  0.08  a~  and  were  highly  variable  with  depth. 

In  contrast  to  the  general  agreement  between  the  estimates  of  absorption  derived  from 
the  models,  ttiere  is  a  large  disparity  between  the  model-derived  values  and  those  measured 
spectrophotometrically  (the  sum  of  the  components),  the  sum  of  the  coaponents  at  441  nm  is 
two-fold  higher  than  the  absorption  values  derived  froa  the  models  for  the  upper  30  m, 
although  agreement  improves  with  increasing  depth,  also,  there  are  many  values  for  520  nm 
absorption  which  fall  below  0.048.|p~^,  the  absorption  ascribed  to  molecular  water,  that  can 
be  attributed  to  raaan  scattering^  . 

Station  11 

Station  11  consisted  of  five  tlae  series  and  exhibited  low  a  values.  Irradiance  data 
were  interpolated  among  the  five  depths.  The  time  series  samples  are  spaced  at  depth 
increaents  that  are  too  large  to  properly  estimate  absorption  coefficients  using  either 
theoretical  technique;  however,  it  is  useful  to  compare  such  estiaates  to  those  derived 
froa  aeasureMnt  of  individual  col^>onent  absorption.  The  absorption  at  441  na  increased 
froa  0.015  a~^  near  the  surface  to  0.045  s?  near  the  fluorescence  aaxiaua.  At  520  na, 
the  total  range  of  a  was  0.040  to  0.044  a  .  Because  the  diffuse  attenuation  coefficients 
were  used  to  calculate  the  irradiances  between  the  holds,  the  three-paraaeter  aodel  and  the 
Preisendorfer/Nobley  aodel  yielded  similar  absorption  coefficients  at  four  of  the  five 
sampling  depths  (Table  1). 

Like  Station  8,  a  froa  the  coa^nent  analysis  was  much  higher  than  that  predicted  froa 
apparent  optical  properties.  The  differences  in  a  were  greater  for  441  na  than  for  520  na 
(Table  1)  and  were  high  near  the  surface  and  at  tKe  fluorescence  maxima.  Only  at  100  a  for 
441  na,  as  the  particle  concentration  was  declining,  was  the  sua  of  component  absorption 
less  than  that  derived  froa  apparent  optical  properties.  Since  dissolved  substances  were 
low  at  Station  11,  the  disparity  in  the  absorbances  could  not  be  attributed  only  to 
difficulty  in  deterainlng  dissolved  substance  absorption.  Even  discounting  dissolved 
substance  absorption,  the  sum  of  particulate  and  "pure  water"  absorption  still  exceeded 
that  derived  froa  the  Preisendorfer/Nobley  inverse  technique.  For  example,  a(520)  derived 
from  the  theoretical  aodel  was  about  10%  below  that  of  "pure  water".  Either  the  aodel  was 
not  accurately  predicting  absorption,  or  absorption  attributed  to  pure  water  was  high,  or 
some  internal  radiance  source,  such  as  Raaan  scattering'^  was  influencing  the  absorption 
derived  from  apparent  optical  properties  (the  theoretical  models). 

4.6.  Errors  in  calculation  of  absorption 

The  differences  between  absorption  values  at  441  na  derived  from  the  theoretical  models 
and  absorption  calculated  as  the  sua  of  the  coaponents  are  indicative  of  errors  due  to 
ove-est last ion  of  absorption  by  dissolved  and  suspended  Mterial.  The  determination  of  a 
is  intrinsically  difficult  and  some  error  is  likely.  In  addition,  overeat laat ion  of  ^ 
particulate  absorption  may  arise  due  to  uncertainty  in  the  path  correction  necessary  for 
the  Nitchell  and  Kiefer  technique.  Neasureaents  from  stations  with  no  detectable  dissolved 
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substance  component  along  with  ancillary  measurements  of  particulate  absorption  using  opal 
glass  cuvettes  have  suggested  the  path  length  factor  may  have  been  underestimated  by  up  to 
40%.  This  would  lead  to  an  underestimate  of  particulate  absorption. 

Absorption  values  derived  from  the  theoretical  models  are  subject  to  uncertainties 
arising  from  internal  radiance  sources  such  as  natural  pigment  fluorescence,  fluorescence 
of  dissolved  substances,  and  Raman  scattering,  and  factors  influencing  in  situ  irradiance 
measurements  such  as  wave  focusing  and  ship  contamination.  The  processing  methodology 
removed  much  of  the  uncertainty  in  irradiances.  Stavn  and  Weidemann^  have  shown  that 
Raman  scattering  can  be  a  very  important  factor  in  explaining  the  disparity  that  can  arise 
between  absorption  values  derived  from  theoretical  models  and  summing  components  of 
absorption. 

4.7.  Particle  Characterisation 


Measurements  of  individual  particle  characteristics  were  made  to  develop  a  data  set  of 
relevant  parameters  from  which  scattering  and  backscattering  coefficients  could  be 
estimated  and  compared  to  the  theoretical  estimates  of  backscatter  derived  from  the 
Preisendorfer/Mobley  model.  The  results  of  the  beam  transmissometry  and  fluorometry,  two 
accepted  in  situ  methods  for  estimating  total  suspended  particle  concentration  and 
phytoplankton  concentrations  respectively,  have  already  been  described.  The  "Coulter 
counter",  or  more  properly  the  resistive  pulse  technique  is  an  accepted  method  for 
estimating  particle  concentration  and  size  (volume)  distributions.  "SAX". XScanning 
electron  microscopy  with  Automated  image  analysis  and  X-ray  spectroscopy)'^  is  an  automated 
method  whose  applicability  to  particle  characterization  is  being  tested.  SAX  estimates 
particle  concentration  and  size  (projected  area)  distribution,  and  additionally  gives  the 
elemental  composition  for  each  particle,  from  which  we  can  deduce  a  refractive  index.  For 
comparison  to  the  resistive  pulse  data,  particle  volumes  were  estimated  from  the  SAX  data 
by  calculating  the  diameter  of  a  circle  with  an  area  equivalent  to  the  measured  particle 
and  then  finding  the  volume  of  an  equivalent  sphere.  Only  the  data  from  Station  11  will  be 
discussed.  Because  the  irradiance  values  were  too  low  to  reliably  estimate  backscatter 
below  100  m,  only  the  top  three  water  sample  depths  (10,  30,  and  75  m)  were  analyzed  for 
particle  characterization. 

The  resistive  pulse  counts  of  the  total  particles  in  the  size  range  of  1  to  100  urn  m 
of  seawater  are  about  an  order  of  magnitude  higher  than  the  SAX  results,  and  the  depth 
profiles  of  particle  counts  also  differ  between  the  two  techniques  (Figure  3).  The 
discrepancy  probably  lies  in  the  lower  portion  of  the  size.Kange.  Particle  concentrations 
increase  logarithmically  with  decreasing  size  in  the  ocean^  ,  and  any  differences  in  the 
way  the  two  techniques  measure  the  smallest  particles  will  tend  to  magnify  the  differences 
in  their  total  counts. 

When  the  particle  numbers  are  converted  to  particle  volume,  the  two  techniques  continue 
to  show  differences.  The  resistive  pulse  technique  detected  particles  as  large  as  128  urn 
at  all  three  depths,  whereas  the  SAX  technique  detected  particles  no  larger  than  40  urn. 

This  is  due  to  the  fact  that  the  seawater  volume  sampled,  300  ml,  contained  insufficient 
numbers  of  large  particles;  however,  the  trend  of  the  differences  between  the  size 
distributions  are  clear:  the  resistive  pulse  data  reveal  a  relatively  flat  distribution 
typical  of  oligotrophic  waters*'  ;  the  SAX  data  show  a  shift  to  fewer  and  larger  particles 
(Figure  5). 

The  added  dimension  of  elemental  analysis  with  SAX  allows  some  categorization  of  the 
particle  population  (Figure  6).  Of  the  total  particle  populations  at  the  three  depths, 
organic  particles  were  the  majority:  52,  44,  and  56%  of  the  particles  at  10,  30,  and  75  m 
respectively.  The  next  most  numerous  category  was  ">95%  silicon",  i.e.  diatoms  and  pieces 
of  diatom  frustules:  16,  27,  and  30%  of  the  population  at  each  depth.  There  was  a  small 
but  persistent  population  of  silicate  particles  at  8,  8,  and  4%  respectively. 

Particles  with  >95%  calcium  made  up  a  small  fraction  of  the  population  at  4,  2,  and  1%,  and 
barite  was  found  in  trace  amounts  at  all  depths.  As  seen  in  Figure  6,  organic  particles 
tend  to  dominate  the  size  distribution  at  all  depths.  At  10  m  silicates  are  more 
significant  than  >95%  Si  particles  in  particles  larger  than  10  urn.  At  30  m  and  75  m, 
organic  and  >95%  Si  distributions  are  about  the  same  up  to  a  diameter  of  about  10  um. 

There  is  a  peak  in  the  >95%  Si  distribution  about  16  um  diazwter  at  both  depths,  and 
organic  particles  apparently  became  larger  with  depth,  peaking  about  20  um  diameter  at  30  m 
and  about  30  um  diameter  at  75  m. 

4.8.  Comparison  of  particle  data  to  backscatter  estimates 

Qualitative  estimates  of  backscatter  based  on  particle  characterisations  differ 
depending  on  which  measure  is  used.  The  beam  transmissometry  data  suggests  that 
backscatter  should  remain  fairly  constant  from  the  surface  to  75  m,  while  the  in  situ 
fluorometry  suggests  that  backscatter  might  increase  slightly  around  75  m  due  to  increased 
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numbers  of  phytoplankton.  Based  on  the  concentration  of  particles  detected  by  the 
resistive  pulse  technique,  backscatter  should  show  a  peak  around  30  m,  but.  based  on  the 
size  distributions,  there  should  be  little  change  in  backscatter  with  depth.  In  contrast, 
the  particle  concentration  detected  by  SAX  suggests  that  backscatter  should  decrease  with 
depth;  however,  the  size  and  composition  data  from  the  SAX  analysis  suggests  that 
backscatter  should  increase  with  depth. 

Backscatter  estimated  from  optical  measurements  using  the  Preisendorfer/  Mobley  model 
shows  an  increase  to  a  maximum  near  75  m  for  441  nm,  while  backscatter  at  520  nm  remained 
fairly  constant  with  depth.  The  relationships  among  the  various  characterizations  of 
suspended  particles  and  the  estimates  of  backscatter  derived  from  the  model  are  not  clear 
at  this  point  due  to  limitations  that  are  inherent  to  the  measurements  and  the  model.  On 
the  one  hand,  the  various  results  of  particle  analyses  reflect  different  aspects  of  the 
population  are  being  measured;  on  the  other  hand,  each  technique  has  its  own  intrinsic 
errors.  For  its  part,  the  Preisendorfer/Mobley  estimates  of  backscatter  at  Station  11 
suffer  from  the  same  limitation  that  was  seen  in  the  absorbance  estimates,  namely  that  the 
model  was  not  intended  for  use  with  only  a  few  widely  spaced  measurement  depths  such  as  the 
5  sampling  depths  spaced  by  nearly  30  m  at  this  station.  In  addition,  backscatter 
estimates  can  be  subject  to  large  uncertainties  brought  about  by  in  situ  radiance  sources 
such  as  Raman  scattering  and  natural  fluorescence. 

5.  CONCLUSIONS 

1.  Models  of  Stavn  and  Preisendorfer/Mobley  yielded  absorption  coefficients  that  were  very 
similar  if  the  distribution  factor  was  calculated  in  the  Preisendorfer/  Mobley  model. 

2.  At  441  and  520  nm  the  sum  of  absorption  due  to  molecular  water,  dissolved  substances, 
and  particulates  was  higher  than  absorption  values  derived  from  theoretical  models. 

3.  Near  the  surface,  the  values  of  absorption  derived  from  the  theoretical  models 
approached  absorption  coefficients  ascribed  to  molecular  water. 

4.  Preliminary  attempts  to  intercompare  particle  characteristics  and  backscatter  estimates 
derived  from  Preisendorfer/Mobley  model  have  been  largely  unsuccessful.  Increases  in 
backscatter  at  441  nm  were  coincident  with  fluorescence  maxima  and  increases  in  organic  and 
siliceous  particle  size,  but  backscatter  at  520  nm  showed  very  little  structure. 

5.  Internal  radiances  sources  should  be  investigated  as  a  possible  source  of  error  in  the 
theoretical  models,  and  with  its  inclusion  differences  in  absorbances  by  the  various 
techniques  may  diminish. 
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ABSORPTION  (N-'  )  AT  520  MC 


Figure  1.  Absorption  coefficients  at  441  nm 
(Figure  lA)  and  520  nm  (Figure  IB)  derived  using 
the  Preisendorfer/Mobley  model  with  calculated 
average  cosines  (-calculated)  and  the  three- 
parameter  model  of  Stavn  (-•■  uK  ). 
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Figure  4.  Size  distribution  of  particles 
by  Resistive  Pulse  for  each  of  the  three 
depths  sas^led.  "DiaMter”  is  the  diaswter 
of  a  sphere  with  a  voIubm  equivalent  to  the 
Measured  particle  voluae;  ”Total  VoIusm*  is 
the  individual  particle  voluae  for  a  given 
particle  diasMter  tiaes  the  total  nuMber  of 
particles  in  that  size  class  per  ■  of 
Seawater  (equivalent  to  PPB  by  voIusm). 


Figure  5.  Size  distribution  of  particles  of 
all  chemical  classes  by  SAX  for  each  of  the 
three  depths  sampled.  "Diameter"  is  the 
diazMter  of  a  sphere  with  a  projected  area 
equivalent  to  the  area  of  the  particle 
measured.  "Volume"  as  in  Figure  4. 
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Figure  6.  Size  distribution  of  particles 
in  the  three  major  chemical  classes. 
"Diameter"  and  "Volume"  as  in  Figure  5. 
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